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We present an interesting connection between Brownian motion and magnetism. We use 
this to determine the distribution of areas enclosed by the path of a particle diffusing on a 
sphere. In addition, we find a bound on the free energy of an arbitrary system of spinless 
bosons in a magnetic field. The work presented here is expected to shed light on polymer 
entanglement, depolarized light scattering, and magnetic behavior of spinless bosons. 



Consider a particle diffusing on a sphere. If the diffusing particle returns to its start- 
ing point at time (3 its path subtends a solid angle Q at the center of the sphere. 
We ask: what is the probability distribution of Q ? This problem comes up if one 
considers a spin- system in a random magnetic field. As is well known, the state 
(up to a phase) of a spin- ^ system can be represented as a point on the Poincare 
sphere. Under the influence of a random Hamiltonian, the state of the system dif- 
fuses on the Poincare sphere. From the work of Berry 1 and others, it is known that 
the system picks up a geometric phase 7 equal to half the solid angle swept out on 
the Poincare sphere. To compute the distribution of geometric phases one is led to 
the question posed above. A closely related problem has already been studied in 
the context of polymer entanglement: 2 given that a Brownian path on the plane is 
closed at time /3, what is the probability that it encloses a given area Al 

In this paper we present a general method of solving these problems by using a con- 
nection between Brownian motion and magnetism. The qualitative idea is to use a 
magnetic field as a "counter," to measure the area enclosed in a Brownian motion. 
We derive a relation between the distribution of areas in a Brownian motion and 
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the partition function of a magnetic system, which can be used to cast light on both 
subjects. Despite its apparent simplicity, this relation does not seem to have been 
noticed or exploited so far. Our main purpose here is to illustrate its usefulness. 
We first discuss the planar problem solved earlier. We then go on to solve the (as 
yet unsolved) problem of diffusion on the sphere. We also exploit the relation to 
learn about the magnetic properties of bosonic systems. Here we recover previously 
known results and arrive at some others. We conclude the paper with a few remarks. 

Let a diffusing particle start from a point on a plane at time r = O. Given that the 
path is closed at time j3 (not necessarily for the first time), what is the conditional 
probability that it encloses a given area Al By "area" we mean the algebraic area, 
including sign. The area enclosed to the left of the diffusing particle counts as posi- 
tive and the area to the right as negative. This problem has been posed and solved 2 
by polymer physicists, since it provides an idealized model for the entanglement of 
polymers. We present a method of solving this simple problem. 

Let {x(r),0 < r < /3, x(0) = x(/3)} be any realization of a closed Brownian path 
on the plane. As is well known, Brownian paths are distributed according to the 
Wiener measure: 3 if f[x(r)} is any functional on paths, the expectation value of / 
is given by 
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In Eq. (1) the functional integrals 4 are over all closed paths (the starting point is 
also integrated over). (We set the diffusion constant equal to half throughout this 
paper.) Let .4[x(r)] be the algebraic area enclosed by the path x{r). Clearly, the 
normalized probability distribution of areas P(A) is given by 

P(A) = (5(A[x(r)] - A)) w . (2) 

The expectation value 4> of any function <p(A) of the area is given by J P(A)<f)(A)d(A). 
As is usual in probability theory we focus on the generating function P(B) of the 
distribution P(A): 

P(B) = = [ P(A)e ieBA dA, (3) 
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which is simply the Fourier transform of P(A). For future convenience we write the 
Fourier transform variable as eB. The distribution P(A) can be recovered from its 
generating function by an inverse Fourier transform. From Eqs. (2) and (3) above 
we find 

P(B) = (e^ BA ) w . (4) 
Notice that BA can be expressed as 

fP - dr 

BA= / A(M)~dT, (5) 
Jo dr 

where A(x) is any vector potential whose curl is a homogeneous magnetic field B. 
Equations (1), (4), and (5) yield 



P (B) = — ^ r , ... z S -■ (6) 
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By inspection of Eq. (6) we arrive at 

P(B) = Z(B)/Z(0), (7) 

where Z(B) is the partition function (Z(B) = Tr{exp[— /3 if (£>)]}) for a quantum 
particle of charge e in a homogeneous magnetic field B at an inverse temperature 
j3. This is the central result of this paper and it relates Brownian motion and 
magnetism. As the reader can easily verify, the relation (7) holds even if there is an 
arbitrary biasing potential. The plane can also be replaced by a sphere or (IR 3 )^, 
the configuration space of N particles in (M 3 ). In the last case, the area of interest 
is the sum of the weighted areas of the projections of the closed Brownian paths 
onto the x— y plane. Now we demonstrate the utility of Eq. (7) by computing 
the distribution of areas for diffusion on a plane. The partition function Z(B) for 
a particle of unit mass in a constant magnetic field, is easily computed from the 
energies E n — (n + \)eB and degeneracy (or the number of states per unit area) 
(eB/2ir) of Landau levels 5 (throughout this paper we set fx = c = 1): 

From (7) we find P(B) = ((3eB /2)[smh.(j3eB /2)}~ 1 . Taking the Fourier transform 
of P(B) by contour integration we get the result P(A) = (7r/2/5)[cosh(7Ty4//5)]~ 2 
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derived in Ref. 2. This provides a check on Eq. (7) and illustrates its use. 

Let us now address the problem posed at the beginning of this paper: what is the 
distribution P(£l) of solid angles enclosed by a diffusing particle on a unit sphere? 
Unlike the planar case, P(£l) is a periodic 6 function with period 4ir. The generating 
function P g of the distribution of solid angles is given by 

/•47T 

P g = dnP(Q)e ia r (8) 
Jo 

with g an integer. P(fl) is expressed in terms of P g by a Fourier series 

oo 
g=-co 

rather than an integral (3). Relation (7) now takes the form 

P 9 = z 9 /z , (10) 

where Z g is the partition function for a particle of charge e on a sphere subject to 
a magnetic field created by a monopole of quantized strength G = g/e (Ref. 7) at 
the center of the sphere. The energy levels of this system are easily computed: 8 

^■ = b'(+i)-s 2 ]A 

where j, the total angular momentum quantum number ranges from \g\ to infinity, 
and the jth level is (2j + l)-fold degenerate. The partition function is consequently 
given by 

oo 

Z 9 =2_^{1] + l)e * . (11) 

i=\g\ 

Combining (9), (10), and (11) and rearranging the summations we arrive at 

i 00 r ii/- 
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where ({I, (3, fi) = exp[-l/2{/?(2Z + 1) + in}]. The function (12) is plotted numeri- 
cally for various values of f3 in Fig. 1. The qualitative nature of these plots is easily 
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understood. For small values of (3 the particle tends to make small excursions and 
its path encloses solid angles close to or 47r and consequently the plots are peaked 
around these two values. As the available time (3 increases, other values of Q are 
also probable and the peaks tend to spread and the curves to flatten out. Finally in 
the limit of (3 — > oo the particle has enough time to enclose all possible solid angles 
with equal probability. These plots give the answer to the question that was raised 
in the beginning of the paper. 

Now we turn to the magnetic properties of spinless bosons. An N particle system 
in three dimensions placed in a homogeneous external magnetic field which is along 
the z direction has the Hamiltonian 
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FIG . 1 . The probability distribution P(Q) of solid angles for closed randon walks lasting a 
time f3. P(Q) is plotted above for four values of (3:0.5, 2, 5, and 10. 

H&,? a ) = £ &- e « A ^ 2 + V (x"), 
(i=i 

where A(x a ) is the vector potential of the external magnetic field. V(x a ) includes 
an arbitrary interaction between the particles as well as an external potential, m a 
and e a are the masses and charges of the particles. {x a , a = 1,2,..., N} are the 
position vectors of the N particles. The configuration space of the system is given 
by Q = (M 3 )^/ ~ where ~ means that we identify points in (M 3 )^ which differ 
by an exchange of identical particles. For simplicity we give the argument for N 
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identical particles with unit mass and charge e a = e. The argument is easily adapted 
to several species of particles of arbitrary charge and mass. 

Now consider a diffusion on Q biased by the potential V(x a ). The Wiener measure 
is now appropriately modified: 



fV[x(r)}f[x(r)]ex V 
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The area whose distribution we are interested in is defined as follows: Let q(r) be a 
closed curve in Q ■ q(r) determines trajectories of particles {x a (r), a = 1,2,..., N} 
in M 3 . The area functional of interest is ^l[g(r)] = Y2 a I A(x a ) ■ dx a . The area 
functional has the following interpretation. If the final positions of the N particles 
are the same as the initial ones (direct processes), v4[g(r)] is simply the sum of the 
areas enclosed by the projection of the particle trajectories on the (x — y) plane. If 
the final positions differ from the initial ones by a permutation (exchange processes), 
the projections of the particle trajectories still define closed curves on the (x — y) 
plane. ^4[g(r)] is defined as the sum of areas enclosed by these closed curves. 

As before we find that P{B), which is the Fourier transform of the distribution 
P(A) = (5(A[x(t)} — A))w(y) of areas, is given by Eq. (7). It is crucial for our 
argument that the particles obey Bose statistics. 9 Since P{A) is a probability dis- 
tribution, P(B) = Z(B)/Z(0) is the Fourier transform of a positive function. This 
places strong restrictions on the partition function Z{B). Let u^i = l,...,n be 
n real numbers. If one defines the n x n matrix D$ = P(ui — uj), the necessary 
and sufficient condition for P(B) to be the Fourier transform of a positive function 
is 10 ' 11 

AW = DetD (n) > for all n. (13) 

This imposes restrictions on the free energy F(B) = — (1//3) In Z(B) of the system 
in the presence of a magnetic field B. 

For the simplest nontrivial case n = 2, the inequality (13) with B = u\ — U2 leads 
to 

Z(B) < Z(0) (14) 

or equivalently, F(B) > F(0). Since the free energy of the system increases in the 
presence of a magnetic field, the material is diamagnetic. This universal diamagnetic 
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behavior of spinless bosons at all temperatures is known in the mathematical physics 
literature. 11 However, our approach may be accessible to a wider community of 
physicists. Our approach relating Brownian motion to magnetism enriches both 
fields and provides each field with intuition derived from the other. For instance, 
the zero-field susceptibility \ — —d 2 F(B)/dB 2 \B=o of the magnetic system is related 
to the variance of the distribution of areas in the diffusion problem: 

X = hnP(B)]"\ B=0 = -j(A-Ay = -jVmA < 0. (15) 

It is curious that the zero-field susceptibility can be interpreted as the variance of 
the distribution of areas. Since the variance cannot be negative, it follows that x, 
the zero-field susceptibility cannot be positive and so these systems are diamagnetic. 

(3) 

Next consider the case n = 3. The 3x3 matrix D)/ will then be a function of 
u = iii — U2 and v = u 2 — u^{u\ — M3, being expressible in terms of u and v). If we 
set u = (i.e., set U\ = u 2 = 0,w 3 = —v), we find that A^ 3 '(u, v)\ u=0 = 0. It then 
follows from the inequality (13) that A^ 3 ^(m, v) has a minimum at u = for all v. This 
implies that d 2 A^/du 2 \ u=0tV=B > 0. Defining the function U{B) = P' 2 [l - P 2 }' 1 , 
where the prime means derivative with respect to the magnetic field, we find 



U(B)<U(0), (16) 

As can be seen by taking the l imit B -> 0, U(0) = -P"(0) = -Px(0). We de- 
fine a critical field B c = n/[2y/ — /3x(0)]. The inequality (16) implies a bound on 
the partition function. Notice that P' lies in a cone defined by the lines of slope 



-(tt/2B c )J(1 - P 2 ) and (tt/2£ c W(1 - P 2 ). It follows that 

P(B) > cos(irB/2B c ) for \B\ < B c . (17) 

The diamagnetic inequality due to Simon and Nelson 11 gives an upper bound on the 
partition function Z(B) of a system of spinless bosons. The new inequality stated in 
(17) gives us a lower bound on Z(B) (see Fig. 2) (or equivalently, an upper bound 
on the free energy). 
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Fig. 2. The region forbidden by the bounds [inequalities (14) and (18) on the partition 
function. These bounds are shown as solid lines. The dotted curve is the partitition for a 
charged simple harmonic oscillator in an external magnetic field. Notice that the dotted 

curve lies outside the forbidden region. 



As an explicit check on this new bound on the free energy we considered a sim- 
ple system-a charged particle in a magnetic field subject to a harmonic oscillator 
potential. The calculated partition function of this system is close to, but above 
the lower bound set by (17). Needless to say, our bound is derived for an arbitrary 
interacting system of spinless bosons. The new bound presented here along with 
the earlier (14) diamagnetic inequality 11 places strong restrictions on the partition 
function of a bosonic system in the presence of a magnetic field. We find a curious 
and immediate consequence of these restrictions: if the zero-field susceptibility of 
the system vanishes, then Eqs. (14) and (17) imply that Z(B) = Z(0), i.e., the 
system is nonmagnetic at all fields. 

The key result of this paper is a connection between two apparently distinct classes 
of problems-Brownian motion and magnetism. This allows us to compute the dis- 
tribution of solid angles enclosed in Brownian motion on a sphere. As mentioned 
earlier, this problem comes up when computing the distribution of Berry phases in 
a random magnetic field. A more classical context is depolarized light scattering. 
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As is well known, a light ray following a space curve picks up a geometric phase, 12 ' 13 
equal to the solid angle swept out by the direction vector. If a light ray inelastically 
scatters off a random medium, its direction vector does a random walk on the unit 
sphere of directions. The distribution P(Q) computed here is relevant to the extent 
of depolarization in such an experiment. 14 

In the domain of magnetism we find an independent way of arriving at the diamag- 
netic inequality 11 which states that the free energy of a system of spinless bosons 
always increases in the presence of a magnetic field. Spinless charged bosonic sys- 
tems occur in the context of superconductors (which are perfect diamagnets) and 
neutron stars. 15 We believe that the community of physicists working in these areas 
may not be aware of the general results available in the mathematical literature. 
For instance, the diamagnetism of bosons may be relevant 16 to the interpretation of 
recent experiments 17 on high-T c superconductivity. 

Throughout this paper we have only discussed homogeneous magnetic fields. It 
is easy to generalize our discussion to take into account arbitrary inhomogeneous 
fields: all one does is consider the distribution of weighted areas. An obvious appli- 
cation of this is the computation of the probability of entanglement of a polymer 
with a background lattice of polymers. We expect the new method outlined here to 
shed light on open problems in polymer entanglement involving more complicated 
configurations of polymers than the simplest one solved so far. One can also use 
the relation (7) to compute the distribution of winding numbers in diffusion in a 
multiply connected space. 

It is a pleasure to thank N. Kumar for bringing up the problem of diffusion on a 
sphere and several discussions on this work; Barry Simon for his help in finding Ref. 
11; Diptiman Sen for discussions and for giving us Ref. 8, and R. Nityananda for 
discussions and drawing our attention to Ref. 10. 
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